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ABSTRACT 

We present a detailed study of the dust energy balance in the Sombrero galaxy Ml 04. From 
a full radiative transfer analysis, including scattering, absorption and thermal re-emission, we 
construct models that can reproduce images at optical/near-infrared wavelengths, the observed 
stellar SED and the minor axis extinction profiles in the V and Rc band. A standard model, 
that contains only an old stellar population to heat the dust, underestimates the observations of 
dust emission at infrared wavelengths by a factor of ~ 3. Supplementing this basic model with 
a young stellar component of low star formation activity in both the inner disk (SFR ~ 0.21 
M0 yr~^) and dust ring (SFR ~ 0.05 yr~^), we are capable of solving the discrepancy 
in the dust energy budget of the Sombrero galaxy at wavelengths shortwards of 100 jim. To 
account for the increased FIR/submm emission beyond 100 /im, we propose a additional dust 
component distributed in quiescent clumps. This model with a clumpy dust structure predicts 
three-quarters of the total dust content (~ 2.8 x 10^ M©) to reside in compact dust clouds 
with no associated embedded sources. Although the assumption of a clumpy dust structure 
in the Sombrero galaxy is supported by high-resolution optical data, we cannot rule out the 
possibility that dust grains with a higher dust emissivity account for part of the discrepancy in 
the energy budget at submm wavelengths. 

Key words: radiative transfer - dust, extinction - galaxies: ISM - infrared: galaxies - galax- 
ies: M104 



1 INTRODUCTION 

Although the dust content of a galaxy is often negligible compared 
to the other components in the ISM, dust grains play a prominent 
role in several astrophysical processes in a galaxy. On the one hand, 
they regulate the cooling of gas in the ISM through photoelectric 
heating and inelastic interactions with gas particles. This process 
of cooling allows the gas in a galaxy to eventually condense into 
stars. On the other hand, the extinction of stellar photons by dust 
particles distorts our view on the stars populating a galaxy. A pro- 
found knowledge of the dust properties is thus of great importance 
to understand galaxy evolution processes such as star formation 
and determine the specific amount of stellar photons reprocessed 
by dust particles. 

Optical studies reveal the extinguishing effect of dust grains 
on stellar light at different wavelengths. Radiative transfer codes 
have been used to model the dust distribution in spiral galaxies 
from optical data (Kylafis & Bahcall 1987; Xilouris et al. 1997, 
1998, 1999; Alton et al. 2004; Bianchi 2007). These optical dust 
studies all show that the scale height of the dust disk is about half 
of the scale height of the stellar disk, while the scale length of this 
dust disk is often more extended than that of the stellar distribution. 
These thinner dust disks are optically thin (central optical depth < 
1) at optical wavelengths up to an inclination of about 60°. 



An alternative way of studying the dust characteristics in a 
galaxy is to probe the dust emission at FIR/submm wavelengths. 
With the recent advent of space (Herschel; Pilbratt et al. 2010) and 
ground based (ALMA) facilities in this particular wavelength do- 
main, the dust in nearby as well as distant galaxies can be traced 
with a significantly improved sensitivity and spatial resolution. 
Tracing the reprocessed stellar light at FIR/submm wavelengths, 
one would expect that the thermal emission of dust exactly com- 
pensates for the amount of photon energy absorbed at UV/optical 
wavelengths. 

Ideally, a complementary study of the dust component from 
optical/UV images and its thermal re-emission at FIR/submm 
wavelengths allows to put strong constraints on the properties of 
dust in a galaxy. Up to now, the energy budget has already been 
probed in several edge-on spirals (Popescu et al. 2000; Popescu et 
al. 2011; Misiriotis et al. 2001; Popescu & Tuffs 2002; Alton et 
al. 2004; Dasyra et al. 2005; Baes et al. 2010). The edge-on view 
of these objects offers the opportunity to trace dust out to large 
radii and, hereby, make an accurate estimate of the scale length and 
height of the dust, which can hardly be distinguished in the face- 
on case. Unfortunately, the high inclination of the disk also implies 
that most of the details of the spiral arms in the stellar disk are 
smeared out and most details of the dust distribution vanish when 
integrating along the line-of- sight. 
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In this paper, we present a detailed radiative transfer analy- 
sis of the Sombrero galaxy (Ml 04). At a distance of only 9.2 Mpc 
(which is the average of measurements from Ford et al. 1996 and 
Ajhar et al. 1997), it is one of the most nearby early-type spirals. 
The Sombrero galaxy has been studied in detail, with a particular 
focus on its supermassive black hole (SMBH) (Kormendy et al. 
1996), rich globular cluster population (e.g., Spitler et al. 2008; 
Cockcroft et al. 2009; Harris et al. 2010; Burkert & Tremaine 
2010), stellar population (e.g., Burkhead 1986; Ford et al. 1996; 
Mould & Spitler 2010) and discrete X-ray sources (e.g., Li et al. 
2007, 2010). Its proximity implies that also the distribution and 
characteristics of dust can be examined in much detail. Its incli- 
nation of ^ 84° combines the advantages of these edge-on spi- 
rals with an insight of the inner regions. This allows to put strong 
constraints on the distribution of stars and dust in the disk. While 
IRAS observations had to settle for a global SED study in the in- 
frared/submm waveband due to its limited angular resolving power 
(Rice et al. 1988; Schmitt et al. 1997), Bendo et al. (2006) could 
separate the emission from different components contributing to 
the dust emission at infrared wavelengths. From their image mod- 
eling procedure on the highest resolution Spitzer/MIPS data (24 
and 70 /xm), they could separate the emission from the nucleus, 
the bulge, the inner disk and the dust ring contributing to the dif- 
ferent infrared wavebands. In the submillimeter band, Vlahakis et 
al. (2008) could disentangle the AGN from the dust ring compo- 
nent in the LABOCA 870 /xm and MAMBO-2 1.2 mm images. 
This wealth of available FIR observations (IRAS: Rice et al. 1988; 
MIPS: Bendo et al. 2006; LABOCA & MAMBO: Vlahakis et al. 
2008) and the very regular and symmetric shape of the dust lane in 
Ml 04 permit a spatially resolved analysis of the dust energy bal- 
ance in this nearby galaxy. 

In the past, the Sombrero dust ring has been probed in many 
related studies. Matsumura & Seki (1989) derived a radial distribu- 
tion of dust in the disk from observations of the polarization along 
the major axis. Knapen et al. (1991) were the first to study the ex- 
tinction law in Ml 04. They found a good agreement with the Galac- 
tic extinction law, except for the B and V band, which required the 
assumption of embedded stars in the dust lane or additional fore- 
ground fight. Emsellem (1995) used fiigfi resolution HRCAM 
V, Rc and Ic data to construct a model for botfi tfie stellar and dust 
component in tfie Sombrero ring. From discrepancies between fiis 
model and observed extinction profiles, fie concluded tfiat scatter- 
ing is an important extinction process and neglecting its effect in 
fiigfily-inclined dust disks may cause a substantial underestimation 
of tfie true dust content in sucfi a galaxy. 

In tfiis paper, we will extend tfie analysis in Emsellem (1995) 
to a full radiative transfer study, including absorption, scattering 
and tfiermal dust re-emission. Continuing on tfie image model- 
ing at infrared wavelengtfis from Bendo et al. (2006), tfiis multi- 
wavelengtfi analysis aims to model tfie cfiaracteristics and distri- 
bution of dust in Ml 04 accounting for botfi tfie extinction and 
dust emission properties. We use tfie SKIRT radiative transfer code 
(Baes et al. 2003, 201 1) to construct a 3D model for tfie stellar and 
dust component in tfie Sombrero galaxy. Subsequently, we compare 
tfie FIR/submm emission predicted by our radiative transfer model 
to tfie observed FIR/submm images and fluxes. Tfiis allows an un- 
precedented analysis of tfie dust energy balance in tfie Sombrero 
galaxy. 

Tfie strategy for constructing tfie best fitting model to tfie op- 
tical/NIR data is explained in §2. In §3, we present tfie results from 
our radiative transfer modeling and discuss tfie dust energy bal- 
ance for tfie various models. §4 gives an overview of previous re- 



Table 1. Observed fluxes for the Sombrero galaxy 



Filter 


A (/im) 


Fu (Jy) 


refi 


FUV 


0.15 


0.0056 ± 0.0008 


1 


NUV 


0.28 


0.0177 ± 0.0025 


1 


B 


0.45 


2.25 ± 0.23 


1 


V 


0.55 


2.76 ± 0.28 


1 


R 


0.66 


3.41 ± 0.34 


1 


I 


0.81 


4.30 ± 0.43 


1 


J 


1.25 


8.06 ±0.81 


1 


H 


1.65 


9.19 ±0.927 


^ 




2.17 






IRAC 3.6 


3.6 


3.94 ± 0.53 




IRAC 4.5 


4.5 


2.31 ±0.32 




IRAC 5.8 


5.8 


1.75 ± 0.22 




IRAC 8.0 


8.0 


1.30 ±0.16 




IRAS 12 


12 


0.74 ±0.19 


2 


MIPS 24 


24 


0.65 ± 0.07 


3 


IRAS 25 


25 


0.50 ±0.13 


2 


IRAS 60 


60 


4.26 ± 1.07 


2 


MIPS 70 


70 


6.7 ± 1.3 


3 


IRAS 100 


100 


22.86 ± 5.72 


2 


MIPS 160 


160 


35.1 ±7.0 


3 


LABOCA 870 


870 


0.924 ± 0.092 


4 


MAMBO 1200 


1200 


0.442 ± 0.044 


4 



^References: Dale et al. 2007 (1), Rice et al. 1988 (2), Bendo et al. 2006 
(3), Vlahakis et al. 2008 (4) 

suits from energy balance studies in edge-on spirals and makes a 
comparison witfi tfie main conclusions for tfie energy budget in tfie 
Sombrero galaxy. Our final conclusions are summarized in §5. 



2 MODELING PROCEDURE 

2.1 Data 

In order to recover a reliable representation of tfie distribution and 
properties of botfi stars and dust in tfie Sombrero ring, we will fit 
a 3D model to tfie stellar SED, tfie images and minor axis extinc- 
tion profiles in botfi tfie V and Rc band. Tfie stellar emission in 
tfie SED is constrained by tfie GALEX FUV and NUV, optical 
BVRI, 2MASS JHK and IRAC 3.6 and 4.5 /xm fluxes, wfiicfi 
were taken from tfie Dale et al. (2007) broadband spectral atlas 
(see Table 1 for an overview). IRAC 5.8 and 8.0 /xm measurements 
are also available, but were not used in tfie fitting procedure since 
tfie emission at tfiose wavelengtfis is not only stellar in origin. Also 
tfie GALEX FUV and NUV data were omitted from tfie fitting 
procedure, since tfiey critically depend on tfie amount of extinction 
and, tfierefore, tfie specific amount of young stars in a galaxy (see 
Section 3.2 later on). Tfie V band images were obtained from tfie 
arcfiive of tfie Spitzer Infrared Nearby Galaxies Survey (Kennicutt 
et al. 2003), wfiile tfie Spitzer 24 and 160 /xm images could be re- 
trieved from tfie Spitzer Local Volume Legacy survey (Dale et al. 
2009). 

2.2 Radiative transfer: SKIRT 

SKIRT is a 3D Monte Carlo radiative transfer code, wfiicfi was ini- 
tially developed to investigate tfie effects of dust extinction on tfie 
pfiotometry and kinematics of galaxies (Baes & Dejongfie 2002; 
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Baes et al. 2003). Over the years, the radiative transfer code evolved 
into a flexible tool that can model the dust extinction, including 
both absorption and scattering, and the thermal re-emission of dust, 
under the assumption of local thermal equilibrium (LTE) (Baes et 
al. 2005 a,b). This LTE version of SKIRT is capable of modeling 
the dust properties in a suite of different environments: circumstel- 
lar disks (Vidal & Baes 2007; Vidal et al. 2011), clumpy tori around 
active galactic nuclei (Stalevski et al. 201 1) and a variety of galaxy 
types (Baes et al. 2010; De Looze et al. 2010; Gadotti et al. 2010; 
Gomez et al. 2010). Recently, the code was adapted to include the 
emission from very small grains and polycyclic aromatic hydrocar- 
bon molecules (Baes et al. 2011). 



2.3 Model specifications 

2.3.1 Stellar distribution 

A model for the emission at optical wavelengths in Ml 04 needs to 
account for the contribution from the bulge, disk, nucleus and the 
halo. Considering that only minor departures from axial symmetry 
have been perceived along the disk and in the central region due 
to dust patches, resembling dust lanes in a stellar bar (Emsellem 
& Ferruit 2000), we assume an axi- symmetric geometry for the 
stellar component in the Sombrero galaxy. We adopt the 15 Multi- 
Gaussian Expansion (MGE) components, both in the V and Rc 
band, reported in Emsellem (1995). Only the global intensity of the 
stellar component is scaled linearly until our models achieve the 
best fit with the stellar SED using a minimization procedure. 

The intrinsic SED of the stars is parametrized by a Maras- 
ton (1998, 2005) single stellar population with an age of 10 Gyr 
and a close to solar metallicity Z = 0.016, based on the chemo- 
evolutionary analysis in Vazdekis et al. (1997). This value for the 
metallicity is somewhat higher than the estimate [Fe/H] ^ —0.5 
from Mould & Spitler (2010), who determined the modal metallic- 
ity value for a field 6.8 effective radii from the center of M104. Con- 
sidering that this field is mainly populated by halo stars and, there- 
fore, constitutes mainly of old metal-poor objects, we consider the 
metallicity estimate from Mould & Spitler (2010) within the same 
order of magnitude of the nearly solar metallicity (Z ^ 0.016) 
adopted in this work. 

Besides the main stellar body of the Sombrero galaxy, we also 
added a nuclear source to the model. At optical wavelengths, the 
nucleus of Ml 04 was modeled as a point source in our radiative 
transfer model. Also at IR/submm wavelengths, a significant por- 
tion of the emission arises from the nucleus. Especially the origin 
for the high 850 /xm compared to the 160 /am flux in the nucleus re- 
mains an unsettled issue (Bendo et al. 2006; Vlahakis et al. 2008). 
In order to account for this peculiar emission from the nucleus at 
submm wavelengths, we also include the nucleus of the Sombrero 
galaxy as a point source emitting brightly in the IR. The shape of 
the SED was obtained through an interpolation between the fluxes 
from the nucleus at IR/submm wavelengths (Bendo et al. 2006; 
Vlahakis et al. 2008). 

2.3.2 Dust distribution 

Our dust model will need a representative parametrization of the 
dust distribution in both the ring and inner disk of Ml 04. Similar to 
the stellar component, the very regular and symmetric dust lane in 
the Sombrero galaxy justifies our choice of an axi-symmetric dust 
geometry for the dust ring. Detailed analyses of the extinction pro- 
files along the vertical axis showed that the dust in the dust lane is 



distributed in a ring with three peaks, which average out to a dou- 
ble ring structure along vertical profiles further from the galaxy's 
center (Knapen et al. 1991; Emsellem 1995). 

We model the dust ring as a linear combination of 4 Gaussian 
functions: 



Pd{R,z) = ^ po,^ exp 



{R-R^) 

2a? 



2 



(1) 

For each Gaussian component, the radial distribution is Gaussian, 
while the density in the vertical direction is assumed to be constant. 
This constant density is justified by the small extent of the dust lane 
in vertical direction (hz ^ 20 pc) (Emsellem 1995). Three of the 
Gaussians are used to fit the peaks in the minor axis extinction pro- 
file, the fourth is used to describe the broader, ring-wide distribu- 
tion of dust. 

The spatial distribution of each individual Gaussian compo- 
nent is determined by three parameters: the radial distance Ri at 
which the maximum value of the Gaussian function is situated, the 
width ai of the Gaussian function and the relative mass contribu- 
tion y9o,i. From the extinction curves, one can already determine 
fairly accurate estimates of the radial distance of the density peaks. 
Therefore, the only unknowns which are allowed to vary in our fit- 
ting procedure are the width ai and central density po,i of each 
Gaussian component. 

Although the bulk of the dust content in the Sombrero galaxy 
is distributed in the symmetric ring, the inner disk also contains an 
amount of dust and young stars. In correspondence to the model 
from Bendo et al. (2006), we represent the inner dust disk as an 
exponential disk with a radial scale length of 2000 pc. The scale 
height is confined to 20 pc, following the height of the dust ring. 
The dust content of the inner disk was determined from a modified 
blackbody fit to the MIPS 70 and 160 /xm fluxes and was set to 4.07 
X 10^ M0. 

Throughout the ring and inner disk of the Sombrero galaxy, we 
assume a single dust mixture, consisting of different dust particles 
with a particular grain size distribution. The abundances, extinction 
and emissivity of the dust mixture are taken from the Draine & Li 
(2007) model, which accurately reproduces the optical properties 
of the dust in our own Galaxy. 

2.4 Model fitting 

For a given set of parameters for the stellar and dust distribution, 
we run simulations with the radiative transfer code SKIRT. SKIRT 
self-consistently calculates the temperature distribution and ther- 
mal emissivity of the dust at every position in the model from the 
local interstellar radiation field (ISRF). The result of each simula- 
tion is the observed SED and a set of images at various wavelengths 
ranging from UV to mm wavebands. In order to find the best fit- 
ting model, we have compared the model results with the observed 
minor axis extinction profiles and images in the V and Rc band 
and the stellar SED. Only the minor axis extinction profiles in the 
V and Rc band were used in the fitting procedure, because only in 
those wavebands the signal-to-noise was sufficient in the entire dust 
lane along the minor axis. Allowing the several parameters to take 
values in a confined parameter space, we could narrow down the 
range of valid parameters and finally obtain the best fitting model 
that simultaneously reproduces both the stellar emission and dust 
extinction (through a least- square fitting procedure). 



/. De Looze et al. 



Table 2. The parameters of the 4 Gaussian components 



Standard model 








Md,ring = 7.5 X 10^ M0 








Gaussian component 




(Ti 


M^/Md 


1 


6100 


480 


0.30 


z 


7Qnn 




U.io 


3 


8750 


180 


0.27 


4 


6850 


1500 


0.25 


Model with embedded SF 








Md,ring=7.8 X 106 Mo 








Gaussian component 








1 


6100 


480 


0.33 


2 


7900 


250 


0.18 


3 


8750 


180 


0.26 


4 


6850 


1500 


0.23 



10"" 



3 MODELING RESULTS 
3.1 Standard model 

In a first set of radiative transfer simulations, we only include the 
contribution from the old stellar population in the Sombrero galaxy. 

The spatial distribution of dust in the Sombrero ring was op- 
timized to fit the minor axis extinction profiles in both the V and 
Rc band. Table 2 summarizes the parameter values for the radial 
distance Ri, the width cr^ of the Gaussian function and the relative 
mass contribution MJM^ for each Gaussian component. The rela- 
tive mass contribution MJM^ is in this case a dimensionless num- 
ber, representative for the relative fraction of the total dust content 
in the Sombrero ring. The total amount of dust in this basic model 
is found to be Md = 7.9 x 10^ M©, of which Md = 7.5 x 10^ 
Mq resides in the dust ring and the other Md = 4.1 x 10^ is 
distributed in the inner disk. Figures 1 and 2 show the ability of 
this basic model to reproduce optical data of the Sombrero galaxy. 
In the upper panel of Figure 1, the observed SED fluxes (see Table 
1) are overlaid on the SED model. The middle and bottom rows 
represent the minor axis extinction profiles in the V and Rc band. 
Figure 2 represent the images in the V , MIPS 24 /xm and MIPS 
160 /im wavebands, where the left and right panels correspond to 
the observed data and model images, respectively. 

For this standard model that only accounts for the stellar con- 
tribution from an old stellar population, we examine the dust en- 
ergy balance in Ml 04. When comparing the FIR/submm emission 
predicted from radiative transfer calculations to the observed quan- 
tities, we find a large discrepancy between the radiative transfer 
model and FIR/submm observations. Our standard model underes- 
timates the true FIR/submm SED emission by a factor of ^ 3. This 
lack of IR emission is also obvious if we look at the 24 and 160 /xm 
SKIRT model images (see Figure 2, two bottom rows). Although 
the basic model could reproduce the optical data, modifications to 
this model are necessary to eliminate the remaining discrepancy in 
the IR emission. 



3.2 Model with embedded star formation 

Relying on the large inconsistency between model and observa- 
tions, we require a substantial modification of the standard model 
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Figure 1. The modeled SED and the minor- axis extinction profiles bands 
for the standard model. The top panel represents the total SED as obtained 
with the SKIRT model (solid black line) overlaid with the observed fluxes 
(see Table 1). The contribution of the dust to the SED is also shown. The 
central and bottom panels show the minor-axis extinction profile in the V 
and Rc, where the solid blue and dashed red lines represent the observed 
and modeled extinction profiles, respectively. 



in the previous section to account for the discrepancies in dust 
emission. A possible explanation for the energy balance problem 
in the Sombrero galaxy is the presence of young stars embedded in 
heavily obscured dust clouds. This additional source of dust heat- 
ing will provide a significant boost to the thermal dust heating, 
while their locally absorbed radiation only has a minor contribution 
to the global extinction of starlight. The presence of an additional 
component of embedded stars in the dust lane was already claimed 
by Knapen et al. (1991) to explain the inconsistencies between his 
modeled and observed B and V band extinction profiles. 
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Figure 2. A comparison of the observed and modeled images for the Sombrero galaxy for the standard model. Top, central and bottom rows represent images 
in the optical V band, the MIPS 24 and MIPS 160 bands. Images in a specific band are scaled accordingly for both observed and modeled data, allowing an 
immediate visual comparative inspection. 



Here, we investigate whether an additional component of 
young stars in our model can account for the missing dust emis- 
sion. To obtain a self-consistent model for the Sombrero galaxy, 
we not only require fitting the thermal emission at FIR/submm 
wavelengths, but also the optical SED, V and Rc band images 
and extinction profiles still need to be reproduced. Considering that 
young stars supply an additional dust heating source, the amount of 
dust will need reconsideration in order to still account for the ob- 
served extinction profiles. In response to an additional component 
of young stars, the extinction profile in the V band might decrease 
(depending on the amount of young stars). Although the emission 
from young stars will augment the number of absorption events 
taken place, the total amount of emitted V band photons grows even 
more, causing a drop in the relative amount of extinction. Neverthe- 
less we expect the changes to be small in the V band and almost 
negligible in the Rc band, since the bulk of emission from these 
young stars peaks at shorter wavelengths (UV). 

We now explore the validity of models constructed from our 
standard model complemented with a young stellar population. 
Therefore, we supplement the standard model from the previous 
section with a star formation component both in the dust ring and 
the inner disk. Vlahakis et al. (2008) concluded that the dominant 
heating sources for the dust ring are located on the inside of the 



ring, based on the more extended dust ring at 870 /xm and 1.2 mm 
compared to the emission in the wavelength range from 5.7 to 70 
/xm. Therefore, we insert ^ 75% of this young stellar population in 
the first component of the dust ring (located closest to the nucleus). 
The remaining 25% resides in the second density peak of the dust 
ring. This distribution more or less corresponds to the ratio of the 
dust masses in the first and second Gaussian component describing 
the dust lane in the Sombrero galaxy. For the distribution of star 
forming regions in Ml 04 in the inner disk, we assume the same 
model describing the spatial distribution of dust in the disk. 

The SED of this young stellar population is parametrized as a 
Starburst SED, generated from data of the Starburst99 library (Lei- 
therer et al. 1999), with the same metallicity (Z = 0.016) as the old 
stellar population (10 Gyr) in the Sombrero galaxy. This starburst 
SSP represents a stellar population with a constant, continuous star 
formation, active during the past 100 Myr. 

To achieve the best fitting model that also accounts for the 
excess in the FIRVsubmm emission, we redo the radiative transfer 
modeling for the standard model, complemented with a SF com- 
ponent in the ring and the disk. Constructing a new library of ra- 
diative transfer models, we determine the parameters of the best 
fitting model through a least- square fitting procedure to the SED 
fluxes and images at 24 and 70 /xm, because the bulk of thermal 
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Table 3. The different SFR tracers and estimates for M104. 
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Figure 3. Same as Figure 1, but for the model with embedded star forma- 
tion. 



dust emission heated by this young stellar population will peak at 
these shorter IR wavelengths. 

The best fitting model complements the standard model with 
a stellar component of modest SF activity in the inner part of the 
dust ring (SFR ^ 0.05 Mq yr~^), while a somewhat higher star 
formation activity resides in the inner disk (SFR ^ 0.21 M© yr~^). 
The dust mass in the dust ring was a little increased (Md = 7.8 x 
10^ M0) to maintain the resemblance between the modeled and 
observed V band extinction curve, while the dust mass in the inner 
disk did not require an additional dust mass. The relative distribu- 
tion of dust in the different components of the ring needed some 
fine-tuning due to the presence of a young stellar component (see 
Table 2). 

Considering that the Sombrero ring contains the majority of 
dust mass, it is rather surprising that the bulk of recent star forma- 
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Kennicutt et al. (2009) 
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Kennicutt et al. (2009) 
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tion activity occurred in the inner disk, close to the nucleus. On 
the other hand, a recent image decomposition analysis based on 
IRAC data supports a more disk-like bulge structure, distinct from 
the prominent stellar halo in Ml 04 (Gadotti et al. 2010). In view 
of this scenario, the presence of a young stellar population in the 
plane of the disk-like bulge already seems more plausible. 

This model with embedded SF is able to fit the 24 and 70 /xm 
observations (see top panel in Figure 3). The young stars also pro- 
vide a heating source for very small grains and PAH molecules. 
Indeed, the typical peaks characteristic for PAH emission show up 
in mid-infrared wavebands in the model SED. This result is con- 
sistent with the low-resolution IRS spectrum taken from the ring 
in Ml 04, which shows the most prominent PAH feature in the 
11.3 /xm band, consistent with star formation activity in the ring 
(Bendo et al. 2006). 

Although the observed FUV and NUV fluxes do not agree 
with this specific amount of stellar activity, this dissimilarity is 
likely a consequence of the highly obscured nature of the embed- 
ding young stars accounting for a significant part of the dust heating 
of grains emitting at 24 and 70 /xm. The total star formation activ- 
ity (SFR ^ 0.26 M0 yr~^) in our model agrees well with the SFR 
estimate obtained from Ha and radio observations of the Sombrero 
galaxy in Kennicutt et al. (2009) and is of the same order of mag- 
nitude as other SFR estimates based on a suite of other diagnostics 
(see Table 3). This makes us confident that the star formation ac- 
tivity in this model was not overestimated. Table 3 summarizes the 
predictions for the SFR from different diagnostics and also men- 
tions the reference to the applied SFR relation. All values of the 
SFR were adjusted to the Kroupa (2001) IMF. 

While this model with embedded star formation can account 
for the dust emission at infrared wavelengths shortwards of 100 
/xm, the SED and MIPS 160 /xm image at the bottom row in Fig- 
ure 4 clearly show that the IR emission at longer wavelengths 
(> 100 /xm) is still significantly underestimated. 



4 DISCUSSION 

The most important result from our analysis is the inability to re- 
produce the IR/submm emission at wavelengths longwards of 100 
/xm. The emission from dust grains in our model with embedded 
star formation underestimates the observed IR/submm emission by 
a factor of 3. A similar energy balance problem has been en- 
countered in other spiral galaxies. In this section we discuss the 
most likely explanation for the discrepancy in the energy balance 
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Figure 4. Same as Figure 2, but for the model with embedded star formation. Also the scaling of the image graphics was adopted from Figure 2. 



of the Sombrero galaxy and compare these results with energy bal- 
ance studies in other objects. 

Several edge-on spirals have been the target of a detailed en- 
ergy balance analysis, among which NGC891 stands out as the 
most extensively studied object. In general, those studies indicate 
that 30% of the UV/optical photons is absorbed and re-emitted in 
the infrared, while the optical accounts for the absorption of only 
10% of the stellar radiation (Popescu et al. 2000; Misiriotis et al. 
2001; Popescu & Tuffs 2002; Alton et al. 2004; Dasyra et al. 2005; 
Baes et al. 2010). This inconsistency is referred to as the 'dust en- 
ergy balance problem' . 

Popescu et al. (2000); Popescu et al. (2011) and Misiriotis et 
al. (2001) tried to explain the discrepancy in the energy budget by 
including a second thin diffuse dust disk, associated with a young 
stellar population. Although this additional dust accounts for the 
missing observed FIR/submm emission, Dasyra et al. (2005) and 
Bianchi (2007) show that this second disk scenario contradicts with 
K band observations. Alternatively, a secondary dust component 
can be added in optically thick clumps distributed throughout the 
dust disk. Neglecting a clumpy dust structure can result in an under- 
estimate of the optical dust mass by 40-50% (Witt & Gordon 1999; 
Alton et al. 2000; Misiriotis & Bianchi 2002). These clumps can 
either have embedded star forming regions or have no associated 
sources. The latter clumps are referred to as 'quiescent'. Bianchi 



(2008) shows that a clumpy structure, encompassing half of the 
dust mass (compared to the diffuse dust component) and embedded 
sources, is able to solve the energy balance problem in NGC 891. 
Also the results in Baes et al. (2010) support the presence of heav- 
ily obscured star formation in clumpy dust clouds, to account for 
the missing infrared radiation in UGC 4754. 

An alternative explanation that does not require the assump- 
tion of an additional dust component was explored in Alton et al. 
(2000); Alton et al. (2004) and Dasyra et al. (2005). They argue 
that an increased emissivity in the submm wavelength range com- 
pared to Galactic values, both in a diffuse component and in denser 
environments, also offers a possible solution to the energy balance 
problem. Assuming an increased emissivity by a factor of 4 and 1.5 
at 850 /im and 1.2mm, respectively, compared to the model predic- 
tion from Draine & Lee (1984) and high latitude observations in the 
Milky Way, allows Alton et al. (2004) to account for the discrep- 
ancy in the energy balance in three edge-on spirals (NGC 4013, 
NGC 4565, NGC 5907). Although many authors have explored dif- 
ferent explanations for this excess in dust emission, due to a lack of 
high-resolution multi- wavelength data one has not come to a gen- 
eral consensus regarding the cause of this contradiction in the dust 
energy balance for several edge-on spirals. 

In the same manner that previous energy balance studies have 
shown to underestimate the IR emission by a factor of 3, the 
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standard model in this work, which was constructed from a fitting 
procedure to the available optical data, lacks a similar amount of 
dust emission in the infrared wavebands. 

Relying on the agreement of our model with the SFR calibra- 
tions and 24 and 70 /xm observations, we argue that an additional 
component of heavily extinguished star formation will not bring a 
solution to the remaining discrepancy. Rather than embedded lo- 
calized sources, it seems that either an additional dust component, 
not associated to star forming regions, or dust grains with a higher 
emissivity at FIR/submm wavelengths, compared to Galactic val- 
ues, need to be invoked to explain the discrepancy in the energy 
budget at longer wavelengths (> 100 /xm). 

The former assumption of a dust component distributed in 
clumps throughout the dust ring is a reasonable assumption. Em- 
sellem (1995) already inferred a clumpy dust structure from his 
absorption Ay map in the southern part of Ml 04 (see Figure 4 in 
his paper). Also, high resolution optical HST data (Christian et al. 
2003) of the Sombrero galaxy hint at a clumpy dust structure. 

In late-type galaxies the discrepancy between model and IR 
observations mainly arises at shorter IR wavelengths (NGC891: 
Bianchi 2008 ; UGC 4754: Baes et al. 2010), where the warm dust 
component dominates the IR emission spectrum, and, hence, the 
presence of embedded star-forming clouds seems to be the most 
likely explanation for the dust energy balance problem in those ob- 
jects. In case of the Sombrero galaxy, an early-type spiral (Sa), the 
majority of compact clumps are not associated to localized sources 
and, therefore, have not been triggered to for stars (yet). These com- 
pact clumps are optically thick to the diffuse UV and optical radia- 
tion field and, therefore, will not have a substantial influence on the 
observed extinction profile, because most of the diffuse UV/optical 
radiation from young stellar objects is locally absorbed. Since these 
quiescent clumps are not associated to any source in a star forming 
region, we could consider them as dust clouds for which the con- 
ditions are not appropriate to initiate star formation. Because these 
clouds might still start forming stars, we expect the density of dust 
in these clumps to be very high. Therefore, these compact struc- 
tures could contain a substantial amount of dust in addition to the 
diffuse dust component, which is distributed more smoothly among 
the dust lane. 

Given the low star formation activity in the Sombrero galaxy 
(SFR '-^ 0.26 M© yr~^) compared to late-type galaxies such as 
NGC 891 (SFR 3 M© yr-\ Popescu et al. 2000; Bianchi 2008) 
and NGC 4565 (SFR - 0.7 M© yr-\ Wu & Cao 2006), it is not 
surprising that Ml 04 overall would contain more quiescent dust 
clumps, with no associated embedded sources. To achieve an esti- 
mate of the total amount of dust residing in dense clumps, we calcu- 
late the dust mass required to explain the dust energy balance prob- 
lem in Ml 04. Therefore, we subtract the emission from the model 
with embedded star formation from the observed fluxes at 100, 160, 
870 and 1200 /xm, respectively. By fitting a single-component mod- 
ified blackbody to the residual fluxes, we obtain a dust mass of Md 
= 2.0 X lO'^ M© at a temperature of Td = 16.2 K. Making the as- 
sumption of an additional dust component distributed in quiescent 
compact clumps throughout the dust lane of the Sombrero galaxy 
can account for the IR emission, also at wavelengths longer than 
100 /im (see Figure 5). 

Bianchi (2008) already showed that only a dust distribution for 
which half of the dust is distributed in clumps can explain the en- 
ergy balance in NGC 891. In the Sombrero galaxy, the total amount 
of dust residing in clumps is estimated to exceed the amount of 
diffusely distributed dust and account for about three-quarters of 
the total dust content (Md = 2.8 x 10^ M©). The dust mass in our 
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Figure 5. Spectral energy distribution (SED) for the Sombrero with an em- 
bedded star formation component, supplemented with a dust component 
Md = 2.0 X 10^ M0, distributed in clumps. The solid line represents the 
modeled SED. The modified black-body, modeling the additional dust com- 
ponent in compact clumps, is shown by the blue dashed-dotted Hne. The 
symbols correspond to observations in the optical and infrared spectrum. 

model is a bit higher than previous estimates of the dust content 
in Ml 04, when scaled to the same distance (9.2 Mpc). Emsellem 
(1995) and Bendo et al. (2006) inferred a dust mass (Md = 7.9 x 
10^ M0) comparable to the dust content in our standard model. 
When taking into account observations at longer wavelengths, a to- 
tal dust content of 1.3 x 10^ M© and 1.5 x 10^ M© with dust tem- 
peratures (Td = 22K and 18.4 K, respectively) was found in Krause 
et al. (2006) and Vlahakis et al. (2008) from radio and submm/mm 
observations, respectively. Recent Herschel observations have esti- 
mated a dust mass of Md ^8.1 x 10^ M© at a temperature of Td 
= 22K (Skibba et al. 201 1). The small dissimilarity with the higher 
dust content in our self-consistent Sombrero dust model is largely 
due to the lower dust temperature for the majority of the dust grains 
in our model (Td = 16.2 K). 



5 CONCLUSIONS 

In this work, we have studied the energy balance of the Sombrero 
galaxy in a full radiative transfer analysis, including scattering, ab- 
sorption and thermal dust re-emission. The Sombrero galaxy bene- 
fits from its 84° inclination in such a way that the inner disk of the 
Sombrero galaxy can be revealed, while most details, and hereby 
information about the dust and stellar properties in the disk, vanish 
when integrating along the line-of- sight in more highly inclined ob- 
jects. Complementing this favorable inclination with the very sym- 
metric dust lane in Ml 04 and the availability of high-resolution FIR 
data, the Sombrero galaxy is considered an ideal study object for a 
detailed study of the energy balance. 

From radiative transfer simulations, we obtain a self- 
consistent model for the Sombrero galaxy, which is able to account 
for the observed stellar SED, the extinction of dust particles in the 
dust lane as well as the images at optical and MIR wavelengths. 
This best fitting model assumes a low-activity star forming stel- 
lar component both in the ring (SFR ^ 0.05 M© yr~^) and inner 
disk (SFR ^ 0.21 M© yr~^) of the Sombrero galaxy to account 
for the excess in emission at 24 and 70 fim. Although this model 
with embedded star formation is in agreement with observations of 
the dust component heated by young stars, it still underestimates 
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the dust emission in wavebands longwards of 100 fim. When as- 
suming the presence of an additional dust component, account- 
ing for three-quarters of the total dust content (Md ^ 2.8 x 10^ 
M0), distributed in quiescent compact clumps throughout the dust 
ring, we can account for the observed excess at FIR/submm wave- 
lengths. The presence of a clumpy dust structure is supported by 
high-resolution HST data, revealing the dumpiness of dust in the 
Sombrero dust ring. In correspondence to the lower star formation 
activity in early- type galaxies compared to late-type objects, the 
dust clumps in the Sombrero galaxy do not seem to be associated 
to recent star formation. 

Although a clumpy dust structure is a reasonable assump- 
tion, we are not capable of completely ruling out the possibility 
that part of the discrepancy might be caused by a dust component 
composed of grains with a higher emissivity in the submm wave 
bands, compared to Galactic properties. Future studies of the en- 
ergy balance in thirty galaxies, for which ongoing Herschel obser- 
vations (HEROES, NHEMESES and FRIEDL programs) soon be- 
come available, will significantly improve the sample statistics and, 
hopefully, shed light on the primary cause (additional dust compo- 
nent or dust with different emissivity properties) for the crisis in the 
energy budget in those objects. 
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